Embedded magnetic phases in (Ga,Fe)N: the key role of growth temperature 
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The local chemistry, structure, and magnetism of (Ga,Fe)N nanocomposites grown by metal 
organic vapor phase epitaxy is studied by high resolution synchrotron x-ray diffraction and 
absorption, transmission electron microscopy, and superconducting quantum interference device 
magnetometry as a function of the growth temperature Tg. Three contributions to the magne- 
tization are identified: i) paramagnetic - originating from dilute and non-interacting Fe^+ ions 
substitutional of Ga, and dominating in layers obtained at the lowest considered Tg (800° C); ii) 
superparamagnetic-like - brought about mainly by ferromagnetic nanocrystals of e— FesN but 
also by 7'-Fe4N and by inclusions of elemental a- and 7-Fe, and prevalent in films obtained 
in the intermediate Tg range; iii) component linear in the magnetic field and associated with 
antiferromagnetic interactions - found to originate from highly nitridated Fej;N (x < 2) phases, 
hke ("-Fe2N, and detected in samples deposited at the highest employed temperature, Tg — 950° C. 
Furthermore, depending on Tg, the Fe-rich nanocrystals segregate towards the sample surface or 
occupy two-dimensional planes perpendicular to the growth direction. 



PACS numbers: 68.55.Nq, 75.50.Pp, 75.75.-c, 81.05.Ea 



I. INTRODUCTION 

The epitaxy of magnetically doped semiconduc- 
tors constitutes a versatile mean of fabricating 
in a self-organized way semiconductor/ferromagnet 
nanocompositeai"— with still widely unexplored but re- 
markable functionalities relevant to spintronics, na- 
noelectronics, photonics, and plasmonics. In these 
nanocomposite materials the presence of robust ferro- 
magnetism correlates with the existence of nanoscale vol- 
umes containing a large density of magnetic cations, that 
is with the formation of condensed magnetic semiconduc- 
tors (CMSs) buried in the host matrix and characterized 
by a high spin ordering temperature."* The aggregation of 
CMSs and, therefore, the ferromagnetism of the resulting 
composite system shows a dramatic dependence on the 
growth conditions and co-doping with shallow impurities. 

In particular, the understanding and control of CMSs 
in tetrahedrally coordinated semiconductor films con- 
taining transition metals (TMs) - typical examples being 
(Ga,Mn)As^ (Ga,Mn)N,'^- (Ge,Mn),I (Zn,Cr)Te,^ and 
(Ga,Fe)N2r2i - have been lately carefully considered 
with the necessary aid of nanoscale characterization tech- 
niques. Indeed, the control over the CMSs formation as 
a function of the fabrication parameters, and the possi- 
bility to reliably produce on demand CMSs with a pre- 
defined size, structure, and distribution in the semicon- 
ductor host, arc fundamental requisites for the exploita- 
tion of these nanostructures in functional devices. At the 
same time these studies draw us nearer to understand the 



origin of the ferromagnetic-like features - persisting up 
to above room temperature (RT) - found in a number of 
semiconductors and oxides J^ii^ 

Dilute zincblende (Ga,Mn)As grown by molecular 
beam epitaxy (MBE) is known to decompose upon an- 
nealing with the formation of embedded MnAs nanocrys- 
tals coherent with the GaAs matrix,*"* and a striking 
spin-battery effect produced by these CMSs has been al- 
ready proven**^ As an example of the critical role played 
by growth parameters, MBE Gei_a;Mna; grown below 
130°C is seen to promote the self-assembling of coherent 
magnetic Mn-rich nanocolumns whilst a higher growth 
temperature leads to the formation of hexagonal GesMna 
nanocrystals buried in the Ge host.*^ 

Following theoretical suggestions,— it has recently 
been demonstrated experimentally that it is possible to 
change the charge state of TM ions in a semiconduct- 
ing matrix and, therefore, the aggregation energy by co- 
doping with shallow donors or acceptorsi^^ In particu- 
lar, it has been proven that in the model case of wurtzite 
(wz) (Ga,Fe)N fabricated by metalorganic vapor phase 
epitaxy (MOVPE) the Fermi-level tuning by co-doping 
with Mg (acceptor in GaN) or Si (donor in GaN) is instru- 
mental in controlling the magnetic ions aggregation. *° 

The same system has been thoroughly analyzed at 
the nanoscale by means of advanced electron microscopy 
as well as by synchrotron-based diffraction and absorp- 
tion techniques. The structural characteristics have then 
been related - together with the growth parameters - to 
the magnetic properties of the material as evidenced by 
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superconducting quantum interference device (SQUID) 
magnetometryiSrli It has been concluded that for a con- 
centration of Fe below its optimized solubility limit 
0.4% of the magnetic ions) the dilute system is predomi- 
nantly paramagnetic (PM). For higher concentrations of 
the magnetic ions (Ga,Fe)N shows either chemical- (inter- 
mediate state) or crystallographic-phase separation. ^"^^ 
In the phase-separated layers a ferromagnetic (FM) be- 
havior persisting far above RT is observed, and has been 
related to the presence of either Fe-rich regions coherent 
with the host GaN (in the intermediate state) or of Fcj-N 
nanocrystals in the GaN matrix. These investigations ap- 
pear to elucidate the microscopic origin of the magnetic 
behavior of (Ga,Fe)N reported by other groupsJ^ii^ 

Along the above mentioned lines, in this work we con- 
sider further the MOVPE (Ga,Fe)N material system and 
we reconstruct the phase diagram of the Fe^^N nanocrys- 
tals buried in GaN as a function of the growth temper- 
ature. Synchrotron radiation x-ray diffraction (SXRD), 
extended fine structure x-ray absorption (EXAFS) and 
x-ray absorption near-edge fine structure (XANES), com- 
bined with high-resolution transmission electron mi- 
croscopy (HRTEM) and SQUID magnetometry allow us 
to detect and to identify particular Fe^jN phases in sam- 
ples fabricated at different growth temperatures Tg as 
well as to establish a correlation between the existence 
of the specific phases and the magnetic response of the 
system. Our results imply, in particular, that self assem- 
bled nanocrystals with a high concentration of the mag- 
netic constituent account for ferromagnetic-like features 
persisting up to above RT. These findings for (Ga,Fe)N 
do not support, therefore, the recent suggestions that 
high temperature ferromagnetism of - the closely related 
- oxides is brought about by spin polarization of de- 
fects, whereas the role of magnetic impurities is to bring 
the Fermi energy to an appropriate position in the band 
gapJ^ 

We find that already a 5% variation in the growth tem- 
perature is critical for the onset of new Fe^jN species and 
we can confirm that an increase in the growth tempera- 
ture promotes the aggregation of the magnetic ions, re- 
sulting in an enhanced density of Fe-rich nanocrystals 
in the matrix and in a consequent increase of the fer- 
romagnetic response of the system. Moreover, we ob- 
serve that while in the low-range of growth temperatures 
the Fe-rich nanoobjects tend to segregate close to the 
sample surface, at higher Tg two-dimensional assemblies 
of nanocrystals form in a reproducible way at different 
depths in the layer, an arrangement expected to have 
a potential as template for the self-aggregation of metal- 
lic nanocolumns.^*' The non-uniform distribution of mag- 
netic aggregates over the film volume here revealed, im- 
plies also that the CMS detection may be challenging 
and, in general, requires a careful examination of the 
whole layer, including the surface and interfacial regions. 

The paper is organized as follows: in the next Section 
we give a concise summary of the MOVPE process em- 
ployed to fabricate the (Ga,Fe)N phase-separated sam- 



ples together with a brief description of the character- 
ization techniques. A table with the relevant samples 
and relative parameters completes this part. The central 
results of this work are reported in Section III and are 
presented in two sub-sections discussing respectively: i) 
the detection, identification and structural properties vs. 
Tg of the different Fe^N nanocrystals in phase-separated 
(Ga,Fe)N, with the distribution of the nanocrystals in 
the sample volume, and ii) the magnetic properties of 
the specific families of Fe^^N phases. In Section IV we 
sum up the main conclusions and the prospects of this 
work. 



II. EXPERIMENTAL PROCEDURE 
A. Growth of (Ga,Fe)N 

We summarize here our study by considering a series 
of wurtzite (Ga,Fe)N samples fabricated by MOVPE in 
an AIXTRON 200 RF horizontal reactor. All structures 
have been deposited on c-plane sapphire substrates with 
TMGa (trimethylgallium) , NH3, and FeCp2 (ferrocene) 
as precursors for, respectively, Ga, N and Fe, and with 
H2 as carrier gas. 

The growth process has been carried out according to 
a well established procedure)^! namely: substrate nitri- 
dation, low temperature (540° C) deposition of a GaN 
nucleation layer (NL), annealing of the NL under NII3 
until recrystallization and the growth of a « 1 /im thick 
device-quality GaN buffer at 1030°C. On the GaN buffer, 
Fe-doped GaN overlayers (w 700 nm thick) have been 
deposited at different Tg ranging from 800° C to 950° C, 
with a V/III ratio of 300 [NH3 and TMGa source flow 
of 1500 standard cubic centimeters per minute (seem) 
and 5 seem, respectively] , with an average growth-rate of 
0.21 nm/s, and the flow-rate of the Fe-precursor set at 
300 seem. During the whole growth process the samples 
have been continuously rotated in order to promote the 
deposition homogeneity, while in situ and on line ellip- 
sometry is employed for the real time control over the 
entire fabrication process. 

The considered samples main parameters, including 
the Fe concentration, are displayed in Table ID 



B. Synchrotron x-ray diffraction — experimental 

Coplanar SXRD measurements have been carried out 
at the Rossendorf Beamline BM20 of the European 
Synchrotron Radiation Facility (ESRF) in Grenoble - 
France, using a photon energy of 10.005 keV. The x-ray 
data correspond to the diffracted intensities in recipro- 
cal space along the sample surface normals. The beam- 
line is equipped with a double-crystal Si(lll) monochro- 
mator with two collimating/focusing mirrors (Si and Pt- 
coating) for rejection of higher harmonics, allowing mea- 
surements in an energy range of 6 to 33 keV. The sym- 
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TABLE I: Considered (Ga,Fe)N samples with the correspond- 
ing growth temperature, Fe concentration as evaluated by sec- 
ondary ions mass spectroscopy (SIMS) as well as concentra- 
tion of the dilute paramagnetic Fe'^"'" ions Xp^3+ and a lower 
limit of the concentration of Fe ions xpcj^ contributing to the 
Fe-rich nanocrystals, as obtained from magnetization data. 



Sample 


Tg 


Fe concentration 








°C 


[lO^o cm-^] 


[10" cm- 


[10" cm-3] 


S690 


800 


1 


3.2 


0.1 


S687 


850 


2 


2.9 


1.7 


S680 


850 


2 


2.7 


1.5 


S987 


900 


4 


2.4 


1.6 


S691 


950 


4 


2.9 


3.2 



metric uj = 29 scans are acquired using a heavy-duty 
6-circle Huber difFractometer and the most intense peaks 
are found for 29 up to 40°. 



C. XAFS — experimental and method 

X-ray absorption fine structure (XAFS) measurements 
at the Fe K edge (7112 eV) are carried out at the 
"GILDA" Italian collaborating research group beam- 
line2^ (BM08) of ESRF under the same experimen- 
tal conditions reported in Ref. |Tl|, collecting both the 
XANES and EXAFS spectra, and employing the follow- 
ing method for the analysis. 

A set of model compounds is established: Fe substi- 
tutional of Ga in GaN (FeGa)ri C-Fe2N,23 e-FegN^ 7- 
Fe^N^^ a-Fe^^ and 7-Fei^ For these input structures 
the XANES absorption spectra are calculated using the 
FDMNES codc^"" while the EXAFS scattering expansion 
signals are computed with the feff8.4 code^^ in order 
to de-correlate the structural results to a specific software 
choice. In both cases muffin-tin potentials and the Hedin- 
Lunqvist approximation for their energy-dependent part 
is used, with self-consistent potential calculation for en- 
hancing the accuracy in the determination of the Fermi 
energy {Ep). 

X-ray polarization is taken into account for Feca 
while unpolarized simulations are conducted for the other 
phases assuming a random orientation of the nanocrys- 
tals in the sample. In addition, for XANES the conver- 
gence of the results is tested against the increasing in the 
input cluster size (>150 atoms) and the method is vali- 
dated by experimental values from Fcca and a-Fe. The 
resulting simulated spectra are then convoluted via an 
energy- dependent function as implemented in fdmnes^I 
plus a Gaussian experimental broadening of 1 eV and fit- 
ted to the normalized XANES experimental data in the 
energy range from -20 to 80 eV relative to Ep with a 
linear combination analysis using the Athena graphical 
interface^ to ifeffit.^-^ 

All the possible combinations with a maximum of three 



spectra per fit (maximum of six fit parameters: ampli- 
tude and energy shift) are tested and the best fit is chosen 
on the basis of the statistics, discarding unphysical 
results. Finally, the XANES results are independently 
checked through the quantitative analysis of the EXAFS 
data where the background-subtracted (via the viper 
prograni^) fc^-weighted fine-structure oscillations, x(fc), 
are fitted in the Fourier-transformed space. 



D. High resolution transmission electron 

microscopy — experimental 

The HRTEM studies are performed on cross-sectional 
samples prepared by standard mechanical polishing fol- 
lowed by Ar+ ion milling at 4 kV for about 1 h. Conven- 
tional diffraction contrast images in bright-field imaging 
mode and high-resolution phase contrast pictures were 
obtained from a JEOL 2011 Fast TEM microscope oper- 
ating at 200 kV and capable of an ultimate point-to-point 
resolution of 0.19 nm and allowing to image lattice fringes 
with a 0.14 nm resolution. 

Additionally, energy dispersive x-ray (EDS) analysis 
has been performed via an Oxford Inca EDS equipped 
with a silicon detector to obtain information on the local 
composition. Selected area electron diffraction (SAED) 
and fast Fourier transform (EFT) procedures are em- 
ployed to study scattering orders and d-spacing for re- 
spectively the larger and the smaller nanocrystals. 

E. SQUID magnetometry — experimental 

The magnetic properties have been investigated in a 
Quantum Design MPMS XL 5 SQUID magnetometer be- 
tween 1.85 and 400 K and up to 50 kOe following the 
methodology described previouslyi^ 

The difference between the magnetization values mea- 
sured up to 50 kOe at 1.8 K and 5 K is employed to 
determine the concentration xp^3+ of paramagnetic Fe'^''" 
ions in the layers.— The lower limit of the concentration 
of Fe ions contributing to the Fe-rich nanocrystals, as well 
as an assessment of their Curie temperature is inferred 
from magnetization curves at higher temperatures. 

Finally, measurements of field cooled (EC) and zero 
field cooled (ZFC) magnetization hint to the influence of 
the growth temperature on the size distribution of the 
nanocrystals. 

III. RESULTS 

A. Fe^iN phases vs. Tg in crystallographically 
separated (Ga,Fe)N 

Before entering into the detailed discussion of our stud- 
ies, we would like to point out that the reproducibility 
of the data has been accurately tested and: i) different 
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FIG. 1: (Color online) SXRD spectra for (Ga,Fe)N layers 
deposited at different growth temperatures. Inset: peak at 
35.3° deconvoluted into two components assigned to diffrac- 
tion maxima (111) of e-FesN and (110) of Q-Fe [experiment 
(dotted line) and fit (smooth line)]. 
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FIG. 2: SXRD spectra for a dilute (Ga,Fe)N sample (S690) 
as grown and upon in situ annealing at Ta = 900° C for 1 h, 
indicating that post-growth annealing does not induce - in the 
SXRD sensitivity range - crystallographic decomposition. 



samples grown under the same conditions have been char- 
acterized, 11) all measurements (SXRD, HRTEM, etc..) 
have been repeated in different runs on the same samples 
and we can conclude that both the (Ga,Fe)N structures 
are stable over time and the formation of different phases 
is reproduced when the growth conditions are fidely repli- 
cated. 

In Fig. [T] we report SXRD diffraction spectra for the 
(Ga,Fe)N samples grown at different temperatures, as 
listed in Table I. For the layer S690 fabricated at 800°C 
we have no evidence of secondary phases and only diffrac- 
tion peaks originating from the sapphire substrate and 
from the GaN matrix are revealed, in agreement with 
HRTEM measurements showing no phase separation. 
Moreover, in order to test the stability of the dilute phase, 
we have annealed the samples up to Ta = 900° C and 
in situ SXRD measurements upon annealing do not re- 
veal the onset of any secondary phases, as reported in 
Fig.m in accord with the behavior of dilute Mn in GaN^ 
and in contrast with (Ga,Mn)As where post-growth an- 
nealing is found to promote the precipitation of MnAs 
nanocrystals.^ 

Moving to a Tg of 850°C (S687) different diffrac- 
tion peaks belonging to secondary phases become evi- 
dent. We have previously reported^ that when grow- 
ing (Ga,Fe)N at this temperature, one dominant Fe-rich 
phase is formed, namely wurzite e-FcsN, for which we 
identify two main peaks corresponding to the (002) and 
the (111) reflexes, respectively. A closer inspection of the 
(lll)-related feature and a fit with two gaussian curves 
centered at 35.2° and 35.4°, gives evidence of the presence 
of the (110) reflex from cubic metallic a-Fe. Moreover, 
the broad feature appearing around 38° is associated to 



the (200) reflex of face centered cubic (fee) 7'-Fe4N, that 
crystallizes in an inverse perovskite structurei^ From 
the position of the peak, we can estimate that these 
nanocrystals are strained. 
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FIG. 3: SXRD for a sample (S874) grown at 950° C evidencing 
the aggregation of (200) 7-Fe in the (Ga,Fe)N layer. 

As the growth temperature is increased to 900°C 
(S987) there is no contribution left from the (110) a-Fe 
phase, and the signal from the (111) of e-FeaN is signif- 
icantly quenched, indicating the reduction in either size 
or density of the specific phase. Furthermore, an intense 
peak is seen at 34°, corresponding to the (121) contribu- 
tion from orthorhombic C-Fe2N. This phase crystallizes 
in the a-Pb02-like structure, where the Fe atoms show 
a slightly distorted hexagonal close packing (hep), also 
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found for e-FegN,^ 

The structural resemblance of e-FeaN and the ^-Fe2N 
is remarkable, as the hep arrangement in e-FcsN is nearly 
retained in C-Fe2N.^ This gives a hint of the likely di- 
rect conversion of phase from e-FeaN into C-Fe2N. The 
diffraction peak from (200) 7'-Fe4N is still present at this 
temperature, but its position is slightly shifted to its bulk 
value. A similar behavior is observed for the diffraction 
from (200) e-FegN (002), shifted from 32.78° to 32.9°. 

At a growth temperature of 950°C (S691) the diffrac- 
tion peak of (200) 7'-Fe4N recedes, indicating the decom- 
position of this fee phase at temperatures above 900° C, 
in agreement with the phase diagram for free standing 
Fea;Nf2i reporting cubic 7'-Fe4N as stable at low temper- 
atures. Only the (002) e-FegN- and the (121) C-Fe2N- 
related diffraction peaks are preserved with a constant 
intensity and position with increasing temperature, sug- 
gesting that at high Tg these two phases and their cor- 
responding orientations, are noticeably stable. Further- 
more, in samples grown at this Tg the peak from (200) 
7-Fe is detected around 41.12°, as reported in Fig. [U in 
agreement with the XAFS data discussed later in this 
Section. 
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FIG. 4: (Color online) Average size vs. Tg of nano crystals in 
the different Fe^^N phases, as determined from SXRD. 

Following the procedure employed previously;^ and 
based on the Williamson-Hall formula methodj— we ob- 
tain the approximate average nanocrystals size from the 
full-width at half maximum (FWHM) of the diffraction 
peaks in the radial (w/26') scans. The FWHM of the 
(002) e-FeaN, of the (200) 7'-Fe4N, and of the (121) C- 
Fe2N diffraction peaks are comparable for samples grown 
at different temperatures, indicating that the average size 
of the corresponding nanocrystals is also constant, as 
summarized in Fig. |4l 

The (111) e-FesN signal intensity is seen to change 



abruptly when comparing the results for the sample 
grown at 850°C to those from the layers fabricated at 
higher temperatures. From the FWHM for this particu- 
lar orientation we can estimate that the nanocrystal aver- 
age size adjusts between 16.5 and 12.0 nm in the consid- 
ered temperature range. At high temperatures, the size 
then remains constant up to 950° C. The size of the a-Fe 
nanocrystals can only be estimated for the sample grown 
at 850°C, where the corresponding diffraction peak can 
easily be resolved and suggests an average size of these 
objects larger than that of the other identified phases, as 
confirmed by the HRTEM images reported in Fig. [HI 




1 2 3 4 5 

Distance, R [A] - without phase correction 



7110 7120 7130 7140 7150 7160 7170 7180 7190 
Energy [eV] 

FIG. 5: (Color online) Normalized XANES spectra (main 
plot) and the amplitude of the Fourier-transforms (inset) of 
the fc^-weighted EXAFS in the range from 2.5 to 10.0 A"^ 
for three samples (points) grown at different temperatures, 
with their relative fits (lines) summarized in Table HIl 

The XAFS study on the (Ga,Fe)N samples fabricated 
at different Tg permits to have a structural description 
of the atomic environment around the absorbing species 
from a local point of view, complementary to SXRD. The 
experimental data are reported in Fig. [S] with the rela- 
tive fits obtained by following the method described in 
Sec. Ill CI Qualitatively an evolution with Tg is visible 
and it is quantitatively confirmed by the results summa- 
rized in Table ini 

In particular, from the XANES analysis - sensitive to 
the occupation site symmetry and multiple scattering ef- 
fects - it is possible to infer how the composition of the 
different phases evolves with increasing Tg-. Feca reduces, 
while e-FeaN increases up to 950°C when the precipita- 
tion is in favor of C-Fe2N and 7-Fe. This behavior is con- 
firmed also by the EXAFS spectra given in the inset to 
Fig. O where the first three main peaks in the fit range 
from i?niin to i?max represent, respectively, the average 
Fe-N, Fe-Fe and Fe-Ga coordination. 

In addition, the signal present at longer distances con- 
firms the high crystallinity and permits to include im- 
portant multiple scattering paths in the fits for a better 
identification of the correct phase. In fact, from the Fe- 
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Fe distances it is possible to distinguish Fe^jN (w 2.75 A) 
from pure Fe phases (w 2.57 A), while the distinction 
between a-Fe and 7-Fe is possible with the different mul- 
tiple scattering paths generated from the body centered 
cubic (bcc) and from the fee structure, respectively. 
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FIG. 6: Transmission electron micrographs of different Fe^N 
phases: a) HRTEM image of a e-FesN nanocrystal; (b) the 
corresponding FFT image, revealing that the d-spacing along 
the growth direction is about 0.216 nm. (c) HRTEM image 
on Q-Fe nanocrystal in sample S687; (d) SAED pattern on the 
enclosing area, (e) HRTEM image of a (^-Fe2N nanocrystal; 
(f ) the corresponding FFT image, revealing that the ci-spacing 
along the growth direction is about 0.211 nm. 

The presence of the different Fea;N phases detected 
with SXRD has been confirmed also by HRTEM mea- 
surements on the considered samples, as reported in 
Fig. ini All the HRTEM images presented here have been 
taken along the [1010] zone axis. 

By using the SAED technique for the larger nanocrys- 
tals and a FFT combined with a subsequent reconstruc- 
tion for the smaller objects, we have studied the foreign 
scattering orders and the rf-spacings along the growth 
direction. As shown in Fig. [DJa), the transitional Moire 
fringes indicate that there is a set of planes parallel to 
the GaN (002) ones with a similar d-spacing inside the 
nanocrystal. The corresponding FFT image shown in 
Fig. [6jb) gives an additional diffraction spot close to 
GaN (002), corresponding to a rf-spacing of 0.217 nm, 
matching the (ioo2 of e-FesN. The e-FesN phase is found 
in all the considered samples, with the exception of the 




FIG. 7: (Color online) TEM images: distribution of the Fe- 
rich nanocrystals with increasing growth temperature, (a) Tg 
= 800° C (S690) - dilute (Ga,Fe)N; (b) Tg = 850° C (S687), 
Fe-rich nanocrystals concentrated in proximity of the samples 
interface solely; Tg — 950° C (S691) - Fe-rich nanocrystals 
segregating in proximity of the sample surface (c) and of the 
interface between the GaN buffer and the Fe-doped layer (d) ; 
(e),(f) Tg = 950°C (S876) - to be compared with the TEM 
images of S691 in (c) and (d): reproducibility in the distri- 
bution of the Fe-rich nanocrystals for different samples grown 
at the same Tg. 



one grown at 800°G (S690, dilute). The phase e-FcsN 
has the closest structure to wurtzite GaN and we can as- 
sume that the formation of e-FcaN is, thus, energetically 
favored. '^^ 

The micrograph displayed in Fig. |6l^c) has been ob- 
tained from the layer grown at 850°C and refers to a 
nanocrystal located in the proximity of the sample sur- 
face. The corresponding SAED pattern in Fig. [BKd) re- 
veals that the rf-spacing of the lattice planes overlap- 
ping the GaN matrix has a value of 0.203 nm, match- 
ing the duo of a-Fe. For values of Tg between 900 and 
950°G, nanocrystals like the one represented in Fig. [SJ^e) 
are found. The FFT image shown in Fig. [6jf), reveals 
that the additional d-spacing is 0.211 nm, corresponding 
to the rfi2i of C-Fe2N. 

It should be underlined here that the size of the 
nanocrystals in the HRTEM images is smaller than the 
average value obtained from SXRD. This discrepancy 
originates from the fact that a cross-sectional TEM spec- 
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TABLE II: Quantitative results of the XAFS analysis (best fits). XANES: composition (x) and energy shift relative to Ef {AE) 
for each structure; EXAFS: average distance (R) and Debye- Waller factor (a^) for the first three coordination shells around the 
absorber. For each phase the coordination numbers are kept to the crystallographic ones and rescaled by the relative fractions 
found by XANES and a global amplitude reduction factor, 5*0, of 0.93(5) as found for Feca- Error bars on the last digit are 
reported in parentheses. 
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FIG. 8: A phase diagram of (Ga,Fe)N as a function of the 
growth temperature. 



imen must be rendered considerably thin, in order to 
achieve a sufficient transparency to the electron beam. 
Therefore, a nanocrystal is usually only partly enclosed 
in the investigated area. At the same time, in low mag- 
nification micrographs on thicker volumes the size of the 
objects becomes comparable to the average value deter- 
mined by the SXRD studies. 

Cross-sectional low-magnification TEM measurements 
permit to observe that while at the lower growth tem- 
peratures the Fe-rich nanoobjects tend to segregate close 
to the sample surface as seen in Fig. [7{b), at higher Tg 
two-dimensional assemblies of nanocrystals form in a re- 
producible way - as proven by comparing Fig. [7|c),(d) 
and Fig. [7Je),(f) - and this arrangement is expected to 
be instrumental as template for the self-aggregation of 
metallic nano columns j25 

Summarizing the SXRD, XAFS and HRTEM findings, 
a phase diagram of the Fe-rich phases formed in (Ga,Fe)N 
as a function of the growth temperature is constructed 
and reported in Fig. |51 showing the dominant phases for 
each temperature interval. 

According to Ref. [s^ when the concentration of the 



interstitial-atoms in the s phase is increased by only 0.05 
atoms/100 Fe, a phase transition from e to C occurs. In 
this process the Fe atoms retain their relative positions 
but there is a slight anisotropic distortion of the e lat- 
tice that reduces the symmetry of the (nano) crystal to 
(^-orthorhombic. The hexagonal unit cell parameter ahex 
of e-FeaN splits into the parameters 6orth and Corth in 
C-Fe2N. Moreover, according to the Fe vs. N phase dia- 
gram the orthorhombic phase contains a higher percent- 
age of nitrogen'^'* compared to the hexagonal one, and 
this guides us to conjecture that the higher the growth 
temperature, the more nitrogen is introduced into the 
system. 

Remarkable is the fact that by increasing the growth 
temperature the (002) e-FesN is preserved, while the 
(111) oriented nanocrystals are not detected. A focused 
study would be necessary to clarify the kinetic processes 
taking place between 850°C and 900°C. Moreover, it 
is still to be clarified whether the fact that the e-FesN 
nanocrystals oriented along the growth direction are sta- 
ble, while the ones lying out of the growth plane are not, 
may be related to differences in surface energy. 

The Fe^jN phases found in our (Ga,Fe)N samples are 
listed in Table IIIIl together with their crystallographic 
structure, lattice parameters, d-spacing for the diffracted 
peaks, and magnetic properties. 

Further focused studies are required in order to clarify 
the kinetic mechanisms of segregation and possibly the 
range of parameters that could allow the selectivity of the 
species in different two-dimensional regions of the doped 
layers. 



B. Magnetic properties of Fe^^N phases 

As reported in Table [1111 the different Fe^jN phases we 
identify in the considered samples are expected to show 
specific magnetic responses. The e-FesN phase, predom- 
inant in the samples grown at 850° C, is ferromagnetic 
with a Curie temperature Tq of 575 K42 

The 7'-Fe4N phase, also present though in lesser 
amount in these layers, is FM too, with a Tq of 750 Kj^ 
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TABLE III: Structural and magnetic parameters of the Fe-rich phases found in the considered (Ga,Fe)N samples. 
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For the samples deposited at temperatures above 850°C, 
the dominant and stable phase becomes C-Fe2N. 

The magnetic response of these (Ga,Fe)N layers is 
quite typical for semiconductors containing TM ions 
at concentration above or close to the solubility lim- 
its. Regardless of the prevailing diamagnetic compo- 
nent from the sapphire substrate - that we compensate 
with the procedure detailed elsewhere^ - the field de- 
pendency of magnetization M{H) is characterized pri- 
marily by a dominant paramagnetic contribution at low 
temperatures from diluted substitutional Fe'^+ ions and 
by a superparamagnetic-like component saturating (rela- 
tively) fast and originating from various magnetically or- 
dered nanocrystals with high Fe content. Among them, 
the FM hexagonal e-Fea^^^N have the highest density, 
according to the SXRD and HRTEM studies discussed 
above. Despite the richness of different phases, it is rel- 
atively straightforward to separate these major compo- 
nents and to treat them - to a large extent - qualitatively. 

We begin by noting that the superparamagnetic-like 
component originates primarily from nanocrystals char- 
acterized by a relatively high magnitude of the spin order- 
ing temperature, so that their magnetization MN(r, iJ) 
can be regarded as temperature independent at very low 
temperatures. This means that a temperature depen- 
dence of the magnetization in this range comes from di- 
lute Fe'^"'" ions, whose properties in GaN have been ex- 
tensively investigated previouslyj^i^ 

Accordingly, the concentration of these ions Xp(,3+ can 
be obtained by fitting gfisSxp^a+NoABsiAT,!!) to the 
difference between the experimental values of the magne- 
tization measured at 1.85 and 5.0 K, where ABs{AT, H) 
is the difference of the corresponding paramagnetic 
BriUouin functions ABs{AT,H) = Bs(1.85K,ff) - 
Bs{5K, H). We consider the spin S — 5/2, the corre- 
sponding Lande factor g — 2.0, and treat Xp^3+ as the 
only fitting parameter. 

The values established in this way are listed in Table U 
for the studied samples and they are then employed to 
calculate the paramagnetic contribution at any temper- 
ature according to M = gfj,BSxp^3+ Bs{T, H), which is 
then subtracted from the experimental data to obtain 



the magnitude of the magnetization Mn (T, H) coming 
from nanocrystals. 

For the layers grown at Tg < 900°C, Mn{T,H) satu- 
rates at all investigated temperatures for a magnetic field 
above ~ 10 kOe, as evidenced in Fig.[9lja), pointing to a 
predominantly ferromagnetic order within the nanocrys- 
tals. The values of saturation magnetization M^**' ob- 
tained in this way when plotted vs. temperature as in 
Fig. Hl^b) allow us to assess the corresponding Tc from 
a fitting of the classical Brillouin function to the exper- 
imental points. Furthermore, assuming a value of the 
magnetic moment 2/LtB per Fe, as in e-FcsNj^ we deter- 
mine the concentration of Fe ions xfcn contributing to 
the Fe-rich nanocrystals, as shown in Table I. 

However, for the samples deposited at Tg > 900°C the 
magnitude of M-^{H) saturates only at relatively high 
temperatures, namely around T > 150 K, whereas at 
low temperatures it shows the sizable contribution of a 
slowly saturating component, as shown in Fig.l^^c) where 
magnetization data acquired at 1.85 K for the layer S691 
are reported. 

This new contribution must arise from magnetically 
coupled objects with a spin arrangement other than ferro- 
magnetic. According to the SXRD measurements previ- 
ously discussed and summarized in Fig. [I] the most likely 
candidate is orthorombic ^-Fe2N, antiferromagnetic be- 
low 9 K,^^ or slowly saturating weakly ferromagnetic be- 
low 30 K.^^ In this case, in order to establish AI§f^^ we 
employ the Arrot plot method. The value of M^'^^{T) 
determined in this way is reported in Fig. [Hid), and is 
seen to differ considerably from that of layers grown at 
lower Tg. 

We are able to approximate the experimental values 
of Af^^*(T) with a single Brillouin function only for T > 
150 K (dashed line in Fig. [HUd))- This points to a lower 
value of Tc = 430 K, indicating a shift of the chemical 
composition of £-Fe3_a;N from FesN (x = 0) for Tg < 
900°C to at most Fca.eN (x ^ 0.4) for Tg > 900°C, as Tc 
of e-Fe3_2:N decreases with increasing nitrogen content4^ 

Moreover, the gradually increasing values of M^'^'(T) 
below T < 150 K, marked as the hatched area in 
Fig. ini^d), indicate the presence of even more diluted e- 
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FIG. 9: (Color online) (a) Magnetic field dependence of the nanocrystals magnetization Mn at selected temperatures for sample 
S687 (Tg = 850° C). Each Mn{H) curve has been measured from the maximum positive to the maximum negative field, only 
and the dotted lines obtained by numerical invertion are guides for the eye. The dashed lines represent the saturation level of 
magnetization at each temperature, (b) Bullets - temperature dependence of the saturation magnetization M^''* obtained from 
panel (a). Dashed line - the Brillouin function for the magnetic moment of 2/iB per Fe atom, (c) Three major contributions 
to the total magnetic signal (thick brown solid) for sample S691 (Tg — 950° C): i) paramagnetic from Fe^^ (thick red dashed), 
ii) high-Tc superparamagnetic-like (thin green solid) from the nanocrystals and iii) slowly saturating component (blue short 
dashed). Also here only a half of the full hysteresis loop was measured and the dotted lines obtained by numerical reflection 
are guides for the eye. (d) Bullets - temperature dependence of M^^* for sample S691. Dashed line - the Brillouin function for 
the magnetic moment of 2/xb per atom Fe. The blue dotted line follows the excess of M^=" over the contribution from high Tc 
ferromagnetic nanocrystals. 



Fe3_2:N nanocrystals with x ranging from 0.5 to 1 and 
with a wide spectrum of Tc- Importantly, since e-Fe3_-2;N 
preserves its crystallographic structure and the changes 
of the lattice parameters are minor in the whole range 
< a; < 1, all various £-Fe3_a;N nanocrystals contribute 
to the same diffraction peak in the SXRD spectrum, and 
are detected there as a single compound. 

We note that the presence of either e-FesN or C-Fe2N, 
characterized by a low spin ordering temperature, does 
not hinder the determination of the x-pg3+ values, as both 
compounds have a rather low magnetic moment of O.l/is 
per Fe atom. Accordingly, the resulting variation of their 
magnetization is small comparing to the changes of the 
Fe'^"'' paramagnetic signal at low temperatures. 



The procedure exemplified above allows us to establish 
the lower limit for the Fe concentration that precipitates 
in the form of various Fe-rich nanocrystals (a;Fejv)i 
that is determined by the magnitude of M^'^' at low tem- 
peratures. We again assume 2/iB per Fe atom, as in the 
dominant e-FcaN. These values are collected in Table U 
and plotted as the function of Tg in Fig. [TU] We see, 
that xfcn consistently increases with Tg, and that the 
growth temperature plays a more crucial role than the 
Fe-precursor flow rate^ in establishing the total value of 

Finally, measurements of FC and ZFC magnetiza- 
tion confirm the superparamagnetic-like behavior of 
A/n(T, as reported in Fig.[TTJ As seen in Fig. [TlTa) . 
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FIG. 10: (Color online) Estimated lower limit for the Fe con- 
centration that precipitates in the form of various Fe-rich 
nanocrystals XFe„ as a function of the growth temperature. 



the layer grown at the lowest temperature (S690) shows 
a minimal spread in Tb, having its maximum at Tb < 
100 K, and accordingly a non-zero coercivity is evident 
only at low temperatures. 

In contrast, for most of the studied layers a broad max- 
imum on the ZFC curve, exemplified in Fig. Illf b). in- 
dicates a wide spread of blocking temperatures (Tb) ~ 
reaching RT - and consequently a broad distribution in 
the volume of the nanocrystals. These high values of Tb 
are responsible for the existence of the open hysteresis 
in the M{H) curves seen in Figs. [Sl[a),(c) and thus of a 
non-zero coercivity. This observation again points to the 
growth temperature as to the key factor in the determina- 
tion of the crystallographic structure, size and chemical 
composition of the Fe-rich nanocrystals. 



IV. CONCLUSIONS 

The previous^r— and present studies allow us to draw 
a number of conclusions concerning the incorporation of 
Fe into GaN and about the resulting magnetic proper- 
ties, expected to be generic for a broad class of magneti- 
cally doped semiconductors and oxides. These materials 
show magnetization consisting typically of two compo- 
nents: i) a paramagnetic contribution appearing at low 
temperatures and with characteristics typical for dilute 
magnetic semiconductors containing weakly interacting 
randomly distributed magnetic moments; ii) a puzzling 
ferromagnetic-like component persisting up to above RT 
but with a value of remanence much smaller than the 
magnitude of saturation magnetization. 

According to SQUID and electron paramagnetic res- 
onance^ measurements on (Ga,Fe)N, the concentration 
of Ga substituting the randomly distributed Fe^"*" ions 
increases with the iron precursor flow rate reaching typi- 




100 200 300 400 
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FIG. 11: (Color online) ZFC and FC curves measured in ap- 
plied magnetic field of 200 Oe for samples grown at (a) 800° C 
and (b) 950° C. 



cally the value of 0.1%. Our results imply that the magni- 
tude of the paramagnetic response and, hence, the den- 
sity of dilute Fe cations, does not virtually depend on 
the growth temperature. However, the incorporation of 
Fe can be enlarged by co-doping with Si donors, shifting 
the solubility limit to higher Fe concentrations^^ 

The presence of ferromagnetic-like features can be con- 
sistently interpreted in terms of crystallographic and/or 
chemical phase separations into nanoscale regions con- 
taining a large density of the magnetic constituent. Our 
extensive SQUID, SXRD, TEM, EXAFS, and XANES 
measurements of MOVPE-grown (Ga,Fe)N indicate that 
at the lowest growth temperature Tg = 800° C, a large 
majority of the Fe ions occupy random Ga-substitutional 
positions. However, in films grown at higher tempera- 
tures, 850 < Tg < 950°C, a considerable variety of Fe- 
rich nanocrystals is formed, differing in the Fe to N ratio. 
In samples deposited at the low end of the Tg range, we 
observe mostly e— FcsN precipitates but also inclusions 
of elemental a- and 7-Fe as well as of 7'-Fe4N. In all 
these materials Tc is well above RT, so that the pres- 
ence of the corresponding nanocrystals explains the ro- 
bust superparamagnetic behavior of (Ga,Fe)N grown at 
Tg > 850°C. 
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With the increase of the growth temperature nanocrys- 
tals of (^-Fe2N form and owing to antiferromagnetic in- 
teractions specific to this compound, the magnetization 
acquires a component hnear in the magnetic field. This 
magnetic response has been previously observed and as- 
signed to the Van Vleck paramagnetism of isolated Fe^+ 
ions. In view of the present findings, however, its inter- 
pretation in terms of antiferromagnetically coupled spins 
in nitrogen rich Fe^jN {x < 2) nanocrystals seems more 
grounded. 

The total amount of Fe ions contributing to the for- 
mation of the Fe-rich nanocrystals is found to increase 
with the lowering of the growth rate and/or with the 
raising of the growth temperature. At the same time, 
however, the size of individual nanocrystals appears not 
to vary with the growth parameters. Furthermore, an- 
nealing of (Ga,Fe)N containing only diluted Fe cations 
does not result in a crystallographic phase separation. 
Altogether, our findings indicate that the aggregation of 
Fe ions occurs by nucleation at the growth front and is 
kinetically limited. Moreover, according to the TEM re- 
sults presented here, the spatial distribution of nanocrys- 
tals is highly non-random. They tend to reside in two- 



dimensional planes, particularly at the film surface and 
at the interface between the GaN buffer and the nomi- 
nally Fe-doped layer. 

As a whole, these findings constitute a significant step 
on the way to control the chemistry and local structure 
of semiconductor/ferromagnetic metal nanocomposites. 
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